Statistical linearizations for stochastically quantized fields 
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Abstract 

The statistical linearization method known in nonlinear mechanics and random vibrations theory has been applied to stochas- 
tically quantized fields in finite temperature. It has been shown that even in its simplest form the method yields convenient 
implicit equations for the self-energy, equivalent to the Dyson-Schwinger equations resulting from the summation of infinite 
number of perturbative diagrams. Three examples have been provided: the quantum anharmonic oscillator, the scalar (f> 4 theory 
in three spatial dimension, and the Bose-Hubbard model. The Ramanujan summation has been used to deal with divergent 
integrals and series. 



I. INTRODUCTION 

The method of stochastic quantization invented in the 
end of seventies [T] has immediately grasped interest 
of several members of the field-theoretical community, 
please see [SHI] for reviews. The method has offered a 
version of the quantum theory in which one has to deal 
with c-number fields (instead of operators or path inte- 
grals) at the expense of introduction of another tempo- 
ral dimension and Langevin forces. A very considerable 
amount of work has been applied to develop the the- 
ory and run simulations. In fact, it is the possibility of 
efficient simulations that has strongly stimulated the re- 
search in this method. Initially, it has been worked out 
as a way to obtain Euclidean vacuum expectation val- 
ues of the (products of) quantum fields. Later, it has 
been demonstrated that the method can also work in the 
Minkowski space-time [5J |B]. What is more, it has been 
shown that the stochastic quantization can also be used 
to obtain transition amplitudes between quantum states 
in simple quantum mechanical systems [8] . 

In spite of the initially vigorous development, the inter- 
est in stochastic quantization diminished considerably in 
the nineties. The present authors believe that this hap- 
pened because of the following reasons: (i) the method 
has been worked out with the hope to offer efficient sim- 
ulations of Yang-Mills fields, especially in quantum chro- 
modynamics (QCD). However, the computational QCD 
has had in its disposal very efficient unrelated methods, 
and stochastic quantization has not seemed to offer much 
more advantage from the computational point of view; 
(ii) there have been serious difficulties in the simula- 
tions in Minkowski (rather than Euclidean) space due 
to the instabilities caused by the Langevin forces; (iii) 
it appears that the non-perturbative analytic or semi- 
analytic methods have been somewhat underdeveloped 
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in the stochastic-quantization context - the only works 
on the subject of which are aware are the variational ap- 
proaches of [9] and [10l E] and Hartree-like technique of 
(the above four papers are virtually forgotten); (iv) 
the method has not been sufficiently tested on simpler 
quantum systems - for instance, one would expect that 
before difficult lattice QCD simulations are launched, the 
stochastic quantization should first demonstrate that it 
allows to get, say, correct values of the ground-state en- 
ergy of the helium atom. 

The objective of the paper is to address the point 
(iii) and demonstrate that the stochastic quantization 
can very easily import well-known non-perturbative tech- 
niques from the non-linear mechanics as well as the 
physics of random vibrations. Here, we use the simplest 
of those methods, namely, that of statistical lineariza- 
tion [T5lI17| . please see also [T5HT§] for reviews. It allows 
one to replace a non-linear term in a Langevin equation 
with such a linear term that the expectation value of the 
resulting error is minimized. However, there exist also 
several other linearization criteria. We mention here the 
minimization of the energy error, equality of mean-square 
energies, and equality of mean-square system functions 
(i.e. the deterministic parts of the right-hand sides of 
corresponding Langevin equations). 

Our analyses provided below may be classified as a 
specfic version of Hartree-like approach. The equations 
of motion in the additional time-like variables are lin- 
earized and the linear parameters which enter them are 
found variationally on assuming a Gaussian distribution 
of fields. Since the Hartree approaches have been dis- 
cussed in multitude of papers, there is a natural ques- 
tion whether it is profitable to investigate it again. The 
present authors believe that it is indeed the case be- 
cause the simplest version of statistical linearization pre- 
sented here is only a necessary first step towards devel- 
oping more sophisticated and non-trivial approximation 
schemes based on the so-called higher-order linearization 
as well as building statistically equivalent non-linear solv- 
able models. Less importantly, we have not found any 
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discussion of a Hartree-like approximation in the con- 
text of stochastic quantization except of |12j where it has 
been introduced, so to say, "by hand" rather than with 
the help of a systematic procedure like the one described 
below. 

The main part of the paper is organized as follows. 
To introduce the technique of statistical linearization in 
the stochastic-quantization context, the quantum anhar- 
monic oscillator in non-zero temperature is considered 
in Section 2. Temperature-dependent corrections to its 
frequency are obtained. In Section 3, we consider the 
Langevin equations for </> 4 theory in three spatial dimen- 
sions and finite temperature. We show that the statis- 
tical linearization technique provides us with the correct 
self-consistent equation for the self-energy. Its zeroth- 
temperature limit is analyzed using the Ramanujan sum- 
mation of divergent series. Section 4 contains an anal- 
ysis of the Bose-Hubbard model from the point of view 
of stochastic quantization. The Langevin equations are 
again linearized and Dyson-Schwinger-like equation for 
the self-energy is obtained. Section 5 contains some con- 
cluding remarks. 



II. ANHARMONIC OSCILLATOR 



It has been proved [TJ [5] that the Euclidean, ground- 
state correlation functions (x{ti)x(t2) ...x(t n )} can be ob- 
tained by taking the limit: 



Um s ^ r0 o{x{ti)x{t2)—x{t n )) v , 

where the subscript 77 denotes the averaging over the 
stochastic forces 77. Taking the limit is equivalent to the 
averaging by using a stationary solution to the Fokker- 
Planck equation corresponding to the Langevin equation 

We look for the stationary (s —> 00) solution of Eq.pj) 
subject to periodic boundary condition in the Euclidean 
time r to include finite temperature: 



x(s, 0) = x(s, j3h). 

By "stationary" solution" of ^ we mean the inhomo- 
geneous part of the general solution corresponding to s 
sufficiently large that the initial values of x(s) give van- 
ishing contribution as they are exponentially damped. 

Due to the boundary conditions, it is reasonable to 
expand: 



Let us consider the anharmonic oscillator described by 
the following Euclidean action (confer [2"U]): 



x(s,t) 



x v {s)u v {t), 



Se — 



in 
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dr, (1) 



Ito where r is the Euclidean time, x is the position of 
the oscillator, m - its mass, w denotes the frequency 
of the free oscillations, /? is inverse temperature, and A 
parametrizes the strength of nonlinearity. 

According to the prescription of the stochastic quanti- 
zation, x should satisfy the equation: 



where 



Uu{t) = -4= exp(2niVT/P%) = —. \== cxp(icj ly r). 
y/pn vpn 

The stochastic force can be expanded in the same way: 



V— — 00 



dx 
ds 



SSe 
dx 



+ r](s,r), 



where S/Sx denotes the functional derivative, and s is 
an additional independent variable of the temporal char- 
acter while rj(s, t) is a Gaussian Langevin "force" which 
satisfies the following conditions: 



(r,(s, <)) = 0, (r,(s, tMs', t')) = 2hS(s - s')S(t - t'). 

(2) 

Explicitly, the pseudo-dynamics in "time" s is gov- 
erned by the following non-linear stochastic diffusion 
equation: 



dx d x o t \ /n\ 

— mw x — m—x + r)[s, t). (3) 



ds 



dr 2 



with (Vu(s)Vp(s')) = {2h)8 v ,- p 6(s - s'). 
The amplitudes x v (s) satisfy: 



-m (ujI + u? ) x v - m — — ^ 



ds 



p,a 



Now, we are in position to look for an equivalent linear 
system. It has to take the form: 



dx 
ds 



- + mfllxv = rj v {s). 



(5) 



3! 



The statistical linearization technique in its simplest 
version consists of subtracting the equations ([4] - [5]) , tak- 
ing the square of the resulting error, summing over all v 
and taking the expectation value. This expectation value 
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should be computed with respect to the distribution func- 
tion of x v obtained by solving the full nonlinear system. 
As it is, however, not available, we approximate that ex- 
pectation value by its value obtained from the linearized 
system Thus, we need to minimize the expression: 



,., -p — a ^i/^U I 



X(jX 



i>.a 



The absolute value in the above expression appears 
naturally because x v are, naturally, complex; they must 
satisfy, however, the condition X- v = x* because x itself 
is real. Differentiating with respect to Q v and equating to 
zero the obtained derivative dE/dti^ gives the following 
implicit equation for Q^: 

X 'Ylip.xr( X -p X P X c X p-p-a) 



fif 



(7) 



For Gaussian distribution of x^, the above ratio of mo- 
ments can be simplified [18] : 



x 
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i\pn x dx^ 



XpX(jX^ L — p_ (j) , (8) 



so that Vl^ depends only on the second moments of ir M 

A 



n 2 



2 i 2 



2/3h 



^ ^ (%i/X— u) • 



(9) 



Now, the system ^ can be solved immediately to give 
for s — > oo: 



\Xi/X— vj 



mVl 2 



Let us now write Q 2 — lu 2 + lu 2 + II. Then the "self- 
energy" II satisfies the self-consistent equation: 



n 



A 



1 



2m/3 ^ ul + uj 2 + n 



(10) 



Performing now the frequency summation we find that 
II must satisfy: 



n 



2my/uj 2 + 



= (l + 2n B ( v / ^+n)) , (11) 



where ub(x) is the Bose- Einstein factor (exp(f3hx) — 

In Fig. 1 there is a shaded contour plot of the depen- 
dence of the li/uj 2 on dimensionless coupling constant 
L = hX/(2mui 3 ) and dimensionless inverse temperature 
b = PTilu. 

Let us observe that the dependence of II on A is non- 
analytic even for zeroth temperature. Thus, the result is, 
in a sense, " non-perturbative" . 



Dependence of frequency shift on coupling constant and inverse temperature 




FIG. 1: Dependence of the Tl/u> 2 on dimensionless coupling 
constant L — hX/(2mcu 3 ) and dimensionless inverse tempera- 
ture b = fiTiui. The lighter regions correspond to larger values 
of the 11/ lu 2 



III. STATISTICAL LINEARIZATION OF NON- 
LINEAR SCALAR FIELD THEORY 

We start with the following Euclidean action for the 
(b i model: 



Se = 



tin 



dr \ (PxLe, 



(12) 



where 



C E = ^(drt) 2 + \{V^f + ^^c 2 ^ + |^ 4 . (13) 
The resulting Langevin equation for <f> reads 



d<j) _ 1 d 2 (j> 
ds c 2 dr 2 



V 2 d>- 



■-^4> 3 + r 1 (s,t,T). (14) 



We impose periodic boundary conditions in the Eu- 
clidean time, <j}(s, 0, r) = cf)(s, f3h, r). Also, as it is conve- 
nient to work with discrete variables, we initially assume 
that <f> is also periodic in each spatial variable, that is, 
(j) is defined in the box of the volume Lg with periodic 
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boundary condition imposed in each spatial direction. At 
the end of calculations we shall take the limit Lq — > oo. 
Thus, we expand: 



S,T. 



' ^ = y//3W ^ ^ exp(-iw A T) exp(ikr)0 Ak , 



A k 



(15) 

where uj\ = 2ir\//3h, k = (2ir / L )(m, n,p), and A, to, 
n, p are integers. 

Upon a similar Fourier decomposition of the stochastic 
forces we obtain: 



'Ak 



ds 

where 



N\u + ?7Ak, 
(16) 



^Ak = ^2 E <? !, A 1 ,k 1 0A 2 ,k 2 ^A-A 1 -A 2 ,k-k 1 -k 2 - (17) 
Ai,A 2 ki,k 2 



The corresponding linear system takes the form: 



?Ak 



ds 



(18) 



Using the same prescription as before, i.e., minimiza- 
tion of the expectation value of the error, we find: 



^Ak = -f + k 2 

c 2 



On writing 



2^hL% 



E (^iki< 



Alki 



-Ai ,-kJ 

(19) 



T 



Ak 



n 



and taking expectation values with respect to the 
Gaussian probability density associated with the linear 
system, we obtain: 



n 



EE 



i 



2/3L 3 ^ w 2 /c 2 + k 2 + m 2 c 2 /ft 2 + it 



(20) 



H k M 

or by taking the limit Lq ~ ^ °Oj 



n 



2(2tt) 3 /3 



3ft E 



<Tfc- 



j 2 /c 2 + k 2 + m 2 c 2 /h 2 + n 



(21) 

This is a self-consistent, temperature-dependent equa- 
tion for the self-energy of the self-interacting non-linear 
scalar field. It agrees with that derived in [ST] (please 
see also [22]). Needless to say, the really serious anal- 
ysis starts precisely at this point. It must involve both 



the mass and coupling constant renormalization, and has 
been performed, e.g., in [23l [24] . It has been found that 
the separation of the zero-temperature and temperature- 
dependent terms as well as determination of the thermal 
mass is surprisingly non-trivial, [25] . In view of the fact 
that the matter has been carefully discussed in the above 
papers, we shall not give any detailed analysis of our own. 
We would only like to investigate briefly the consequences 
of application to Eq. (21 1 of the so-called Ramanujan 



technique of summation of divergent series. But before 
doing this we would like to observe that, while there is 
nothing new in Eq. (21), we have obtained it practically 



effortlessly using the simplest, and actually trivial, ver- 
sion of statistical linearization of stochastically quantized 
theory. No diagrammatic analyses or functional differen- 
tiation or integration have been necessary. 

Let us perform the frequency summation as well as 



the angular integration in Eq. (21). The expression for 
IT takes now the form: 



n 



ghc 



k 2 dk 



m 2 c 2 /h 2 + n 



1 + 2n B (c/3hJk 2 + m 2 c 2 /h 2 + IT) ) . (22) 



The integral of the first term on the right-hand side is 
obviously quadratically divergent. 

Let us have a closer look at the limit f3 — > oo. The 
second term (which contains a convergent integral) on the 
right-hand side vanishes in that limit. The self-consistent 
expression for the self-energy takes then the form: 



n = 



ghc 

87T 2 



k 2 dk 



^Jk 2 + (m 2 c 2 /h 2 )+Tl 
ghep 2 



K 2 dn 



(23) 



8^ 2 Jo >/k*+T 

where p 2 = (m 2 c 2 /h 2 ) +IT and we assume that p 2 > 0. 
Euler's substitution \J k 2 + 1 = — k + u yields: 



n 



ghep 2 
32tt 2 J{ 
ghep 2 
= 32tt 2 



{u 2 - l) 2 



udu — 2 



du 
u 



(24) 



In order to deal with the above divergent integrals let 
us invoke the Abel-Plana formula 1251: 



n OO OO ^| n Q 

/ f(x)dx = ^f(n)--f(p)-i 
Jo n=0 1 Jo 



/(it) - f(-it) 



dt. 



(25) 

It makes sense if both the integrals and the series con- 
verge. However, let us make use of the following defini- 
tion. We say that an integral J °° f(x)dx is summable in 
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the Ramanujan sense if the sum of the series Y^^=o /( n ) 
in the Ramanujan sense [27j[28] exists, and the integral 



(f(it) - f(-it))/(e™ - l)dt 



converges in the ordinary sense. Then we write: 



(R) 



7/(0)- 



n=0 



/(it) - /(-it) 



e 27rt - 1 

with obvious meaning of the superscript (R). 

What is called here the "Ramanujan summation" of di- 
vergent series consists in the following. Let ^L>j a(n) be 
a formal (divergent) series, and let R{x) — J2 n >o a ( n + x )- 
Then, if R(x) satisfies the difference equation R(x) — 
R(x + 1) = a(x), the value R(l) gives the sum of the 
series J2n>i a ( n )- The sum is unique if we additionally 

require fulfillment of the condition J. 2 R(x) = and a(n) 
does not grow too fast with n. In [27] an algorithm to 
compute R(x) in terms of the Borel transform of a(n) as 
a well as a useful table of the Ramanujan sums of some 
divergent series is given. 

With the above definitions in mind, we can investigate 
the question of summability of divergent integrals in Eq. 
(24 1 in the Ramanujan sense. 



We find immediately: 



<*> / udu = W 



1 1 

udu-- = --+^ y, 

2 2 ^ 

n=0 



-dt. 



(27) 



As the last integral is equal to 1/12, we have: 



Using now the fact that in the Ramanujan sense the 
sum of the harmonic series is the Euler-Mascheroni con- 
stant 7 while 

WTzW^) dt = ^ {z) + ^ ~ l09{z) > 

where ip(z) is the digamma function |29| (please see 
the formula 6.3.21 of [SH]), we find that, quite trivially, 

dt. 

(30) 



(R) 



0, 



because = — 7- 

As a result of the Ramanujan summation, we get in 
the zeroth-temperature limit: 



n 



64tt 2 



ghc ( m 2 c 2 



n 



3rc 



n 



,hc m 2 c 2 



(31) 



(32) 



64tt 2 h 2 
with g ren =g/(l- (ghc/ Mir 2 )). 

Thus the "renormalization" of the coupling constant 
appears to be finite (provided that ghc ^ 647r 2 ), and the 
same is true about the correction to the mass. This is 
possible only because we have performed a summation 
of two divergent integrals. In the present case, the Ra- 
manujan approach resulted in their trivial elimination. 

We would like to repeat that the formal manipulation 
with divergent series in the spirit of Ramanujan cannot 
be thought of as a substitution of the detailed and careful 
renormalization procedure. Nonetheless, the former is 
perhaps of non-vanishing interest and may be worth of 
some further study. 



(R) 



udu — 



_(R) 



12 



n=0 



(28) 



But, according to [2"7] . 



(R) 



1 — -Bfc+i 
k + 1 : 



where 

(R) 



n = 



are 
5/12. 



the 
and 



(R) 



Bernoulli numbers. 



udu = 0. 



Hence 



What is more, 



(R) 



du 



E; 

n=0 



(< 2 + l)(e 



27Tt 



- 1 

(29) 



)dt. 



IV. BOSE-HUBBARD MODEL 

Let us now consider the one-dimensional Bose- 
Hubbard model. It is of some interest because of the 
potential application of cold Bose gases on lattices as an 
"analog computer" to simulate properties of gauge fields. 

The Bose-Hubbard model describes the 
(thcrmo) dynamics of Bose particles on a lattice under 
the assumption that the particles can jump only between 
neighboring sites, and only the on-site interactions are 
present. Mathematically, it is defined by the following 
(real-time) action: 



S 



ih 
~2 



dou 



dt 



da *m 

dt 



-J I 



fj, c o, rn o. ri 



V 



a* 2 a 2 



(33) 
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where J is the so-called hopping parameter, V is the 
interaction strength, and fic is the chemical potential. 
The index m enumerates the sites on one-dimensional 
lattice, and we assume that size of the lattice is N. 

In the present case, it more convenient to work in real 
time, and perform the analytic continuation to imagi- 
nary time only at the end of calculations. Let us notice 
that the well-known difficulties in the simulations of the 
Langevin equations in the real time do not concern us 
here because the whole procedure is analytic (except of 
the last step). 

Using the prescription by [5] and [BJ, we write 



dam SS , „ da* SS . , . 

tY^+VmM, -^=i— + r,* m (s,t). 



ds Sat 



ds Sa r , 



. (34) 

The complex amplitudes are, however, to be considered 
as independent quantities, so that we effectively complex- 
ify the theory. 

Taking the functional derivatives gives us the system: 

dai dai , 2 

— — = -n—-+iJ{ai + i+ai-i)+ifi c Oii-iVa l a l -eai+r]i, 
as at 

(35) 



da? +da? . T , , . . . , , T7 . . 

-j£ = b^^A a *+i+ a *-l)+ t l jl '> a l -iVa? z ai-ea\ +r ll . 

(36) 

In the above equations e denotes an infinitesimal pos- 
itive constant introduced to guarantee the existence of 
stationary solutions. The stochastic forces satisfy the re- 
lation: 



(rn(s,t)rit,(s' \t') = 2h6 w S(s - s')S(t - t'). (37) 

We impose upon ai and a* periodic boundary condi- 
tions in time t with the period To. We shall also assume 
that the whole lattice is periodic, that is ai = ajq+i for 
some integer N. This allows us to diagonalize the linear 
parts of the above equations by expanding: 



N-l oo 



where (ffc M (s)^, M ,(a')) = ^hk'S^'S(s- s'), and P k 
2nk/N. 

The equivalent linear systems read: 



dfkfj. 
ds 



dft 



k 1 1 



ds 



(40) 



(41) 



where A kfl and A*, are not complex-conjugated quan- 
tities. 

Minimization of the error leads now to the following 
self-consistent equation: 



A k p, = -Jujjfi - 2J cos(P fc ) - fi c - + tj^j (42) 



pis 



or, 



A k p. = -huifj, - 2Jcos(P fe ) - fi c - ie + Il, (43) 
where 



n = 1 . 

iNT LI - hu v - 2Jcos(P p ) - fx c - ie 



(44) 



We now make the rotation in the complex plane to 
obtain the results in imaginary time and write simply 
T = —ifih. Summation over v can be easily performed 
to give the self-energy as a sum over finite number of 
discrete momenta: 



V 



N-l 



n= - E coth((/3/2)(n- 2Jcos(P p ) - Mc )) . (45) 



p=0 



Oil 



v^fE E exp{2mkl/N)exp(-2TTifit/T )f k ^s). 



fc=0 p= — oo 



and analogously for a*, r)i, rj?. 

The complex amplitudes f k[L satisfy the relations: 



d Ik 

' " ' 1 = [ihuifj, + 2iJ cos(P fc ) + i[i c - e] f kft 



ds 



. V 



fpvfqpfk+p-q,fj,+u-p + f/c/(38) 



Let us notice that the zeroth-temperature limit is not 
obvious because we do not know the sign of the expres- 
sion inside the hyperbolic cotangent function. The so- 
lutions to the above self-consistent equation for II pro- 
vide us, upon the analytic continuation from Matsubara's 
to the retarded Green function, with the temperature- 
dependent approximation to the self-energy of a particle 
in the Bose-Hubbard model. While the analysis of Eq. 
( 45 1 is beyond the scope of the present work, we provide 



p,q v,p 



dft 



Fig. 2 which illustrates the general structure of solu- 
tions. In this figure we have plotted curves representing 
the function 



fe/i 



ds 



[ihup + 2iJ cos(Pfc) + i[i c - e] ft 



k ft 



. v 



N-l 



E E fUwh-p+w-u+p + cum v = fix) = ~ e coth (ww* - 2 i cos ( p p) - m )) 



p,q v,p 



p=0 
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FIG. 2: Graphical solution to the self-consistent equation for 
the self-energy in the Bose-Hubbard model for J/V = 0.1, 
H/V = 0.5, N = 10. The solutions (i.e. the values of U/V 
which satisfy ( 45 1 ) are given by the points where the curves 
are crossed by the straight line; (a) fiV = 1, (b) /3V = 20 



(i.e., the right-hand side of (45 1), where v = /3V, 
j = J/V, u = Hc/V as dependent on x = H/V, with 
superimposed straight line y = x. The points at which 
that straight line crosses the curves represent solutions 
of (45 1. The parameters chosen to plot Fig. 2 have been: 
j = 0.1, u = 0.5, v = 1 and 20. 



It is to be noticed that one of the parameters we used 
to plot Fig. 2, namely j = J/V has been taken to be 
0.1. Such a value is quite realistic. However, it obvi- 
ously means that we are in the strong coupling regime, 
for which the statistical linearization may lead to erro- 
neous results. Thus, the above figures are to be under- 
stood merely as an illustration of possible solutions of 
( 45 1 . There is, unfortunately, no general theory which 



could determine whether or not the statistical lineariza- 
tion provides valid approximations in the strong coupling 
regime. 



As is clearly seen from Eq. (|45j) and from the figures, 
the gap equation possesses multiple solutions. A similar 
situation occcurs in both QCD |33] and in the supercon- 
ductivity theory [34] . We plan to discuss the important 
and interesting question of stability of various solutions 
elsewhere. 



V. CONCLUDING REMARKS 

In this paper we have applied the statistical lineariza- 
tion technique to the stochastically quantized field equa- 
tions in Euclidean setting. Three examples have been 
given: the quantum anharmonic oscillator, the real scalar 
<fi 4 theory, and the Bose-Hubbard model. In all three 
cases a self-consistent non-perturbative expression for the 
temperature-dependent self-energy has been provided. 
Only the simplest, Gaussian, version of the linearization 
technique have been employed. As a matter of fact, the 
statistical linearization methods are very rich. At least 
two ways can be used in order to improve the acuracy 
of results. Firstly, one can apply the so-called higher- 
order linearization [TH1 HH [3T]. Within that approach, 
one replaces non-linear terms in the field equations by a 
new dependent variable, which satisfies its own (compli- 
cated) differential equation; the latter is then linearized 
by error minimization. That method is somewhat anal- 
ogous to, but not identical with, the breaking of the 
Schwinger-Dyson (or, analogously, BBGKY) hierarchy 
of the integro-differential equations for the correlation 
functions. Secondly, one can also employ the so-called 
"equivalent non-linearization" [THl HH1 E2] approach in 
which the original non-linear equations are replaced by 
simpler (though still non-linear) equations with known 
statistical properties of solutions. Work is in progress on 
application of both those techniques in the physics of cold 
gases, quantum electrodynamics and quantum optics. 
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